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Abstract We investigate the Josephson effect in the 
graphene nanoribbons of length L smaller than the su- 
perconducting coherence length and an arbitrary width 
W . We find that in contrast to an ordinary supercon- 
ducting quantum point contact (SQPC) the critical su- 
percurrent I c is not quantized for the nanoribbons with 
smooth and armchair edges. For a low concentration 
of the carriers I c decreases monotonically with lower- 
ing W/L and tends to a constant minimum for a nar- 
row nanoribbon with W < L. The minimum I c is zero 
for the smooth edges but eAo/h for the armchair edges. 
At higher concentrations of the carriers this monotonic 
variation acquires a series of peaks. Further analysis of 
the current-phase relation and the Josephson coupling 
strength I c Rn in terms of W/L and the concentration 
of carriers revels significant differences with those of an 
ordinary SQPC. On the other hand for a zigzag nanorib- 
bon we find that, similar to an ordinary SQPC, I c is 
quantized but to the half- integer values (n+l/2)4eAo/H. 



1 Introduction 

Graphene, a single layer graphite, is the two dimensional 
honeycomb lattice of carbon atoms. Recently experimen- 
tal realization of graphene has introduced a new material 
with unique properties [U^GI] which provide the possi- 
bility for designing new carbon based nanodevices. Most 
of the peculiarities come from the electronic structure 
of graphene which is fundamentally different from that 
of a metal or a semiconductor. Graphene has a gapless 
semi-metallic band structure with a linear dispersion re- 
lation of low- lying excitations. This makes the electrons 
in graphene to behave identical to two dimensional mass- 
less Dirac fermions [HG3[5]. The low energy physics of 
electrons is then described by a two dimensional Dirac 
Hamiltonian acting on a four-component spinors in two 
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sublattice spaces of the honeycomb lattice and two in- 
dependent K and K' points in the corresponding recip- 
rocal lattice [TUB]- As a condensed matter counterpart 
of relativistic particles [9], graphene makes possible to 
observe early predictions like Klein paradox and zitter- 
bewegung effect [10]. Also it has attracted an intense 
theoretical and experimental attention to study the ef- 
fects of the relativistic-like dynamics of electrons on the 
different quantum transport phenomena which are al- 
ready known in ordinary conducting systems. 

Already several quantum transport phenomena in- 
cluding the integer quantum Hall effect [2Tl3llTi]. con- 
ductance quantization [12] , quantum shot noise [13] , and 
quantum tunneling 10J have been revisited in graphene 
and found to have anomalous features due to the rela- 
tivistic like dynamic of electrons. Recently the possibility 
of the Josephson coupling of two superconductors by a 
graphene layer was considered by Titov and Beenakker 
[14]. Using Dirac-Bogoliubov-de Gennes formalism [T5] . 
they found that in a ballistic graphene a Josephson cur- 
rent can flow even in the limit of zero concentration of 
the carriers i.e. at the Dirac point. They found a current- 
phase relation I((f) oc cos(0/2)arctanh(sin(0/2)) with 
critical current depending inversely on the junction length 
L. These results for a wide ballistic graphene junction at 
the Dirac point are formally identical to those of an or- 
dinary disordered normal metal [16U17) . 

The graphene Josephson junction has been studied 
experimentally very recently by Heersche et al. |18j , and 
also independently Shailos et al. [19]. They have con- 
firmed the proximity induced superconductivity in the 
graphene samples with superconducting electrods at the 
top of them. Heersche et al. have detected a Josephson 
supercurrent in graphene monolayers at Dirac point with 
a bipolar property due to the electron-hole symmetry 

eg. ' 

While the experiments and theory mentioned above 
had been focused on the wide graphene Josephson junc- 
tions, fabrication of graphene nanoribbons with finite 
widths were reported in a very recent experiment [20j . 
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Already several theoretical works have predicted impor- 
tant properties for graphene nanoribbons resulting from 
the electronic states formed at the edges [2T1I221I23] . So 
one would expect interesting properties for the graphene 
based junctions with different edges which introduce a 
new class of Josephson nanodevices. 

In an ordinary superconducting quantum point con- 
tact (SQPC) shorter than superconducting coherence 
length it was predicted by Beenakker and van Houten 
[24] that the critical (maximum) supercurrent shows a 
step-like variation with width of the contacts of order of 
the Fermi wavelength of electrons. Furusaki et al. [25] 
considering such a system under more general assump- 
tions obtained the same results. The quantization of the 
supercurrent resembles the normal conductance quan- 
tization in a ballistic quantum point contact to 2e 2 /h 
times the number of the transparent quantum channels 
[2"o1l2"T] . The quantization of the supercurrent in units 
of eAo/H (Aq is the superconducting gap in the elec- 
trodes) was confirmed experimentally in semiconducting 
heterostructures [28] , 

In this paper we study the Josephson current in a 
graphene based SQPC with length L smaller than the 
superconducting coherence length £ = hvp/A , and an 
arbitrary width W, with smooth, armchair and zigzag 
edges. Within the formalism of Ref. [15] we find that in 
contrast to an ordinary SQPC 24,25] the critical super- 
current I c in smooth and armchair ribbons with a low 
concentration of the carriers is not quantized but rather 
shows a monotonic decrease by lowering W . For a narrow 
strip W < L the supercurrent reduces to a constant min- 
imum which is eAo / h for armchair edges and vanishingly 
small for smooth edges. We further analyze the product 
I c Rn (Rn being the normal resistance of the SQPC) 
representing the strength of the Josephson coupling, and 
the current-phase relations as a function of W/L. While 
the product I c Rn has a constant value of 2.08Z\o/e for 
a wide contact W/L ^> 1 irrespective of the edges type 
[14] . in narrow junctions the Josephson coupling value 
depends on the type of the edges. By decreasing the 
width in smooth case, the product I c Rn decreases mono- 
tonically and reach a minimum (n/2)Ao/e. However, 
in armchair case it increases and tends to maximum of 
nAo/e for a narrow contact W < L. For graphene SQPC 
the dependence of the current-phase relation on W/L is 
also different from an ordinary SQPC. For sufficiently 
narrow nanoribbons, while the smooth-edge case have a 
sinusoidal current-phase relation like the tunneling con- 
tacts, the armchair-edge junction current-phase relation 
is / = 7 c sin(</>/2) like an ordinary SQPC. 

Far from the Dirac point at high carrier concentra- 
tions, the behavior of I c versus W consists of a similar 
monotonic variation and a series of peaks with distances 
inversely proportional to the chemical potential /i. We 
explain the absence of the supercurrent quantization as 
a result of the significant contribution of the evanescent 
modes in the supercurrent which is a unique property 
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Fig. 1 Schematic of the studied junction structures. Upper 
(lower) panel shows a junction which has zigzag (armchair) 
edges. 



of graphene monolayers. Correspondingly the product 
I c Rn acquires an oscillatory behavior as a function of 
(j,W/Kvf with a period 2ir. The oscillations for smooth 
and armchair edges have a it phase shift due to the dif- 
ference in the zero-state longitudinal momenta of two 
cases. 

For a junction with zigzag edges, because of the valley 
filtering characteristic [23] of the wave functions, situa- 
tion is drastically different from other two edges. We find 
that a zigzag nanoribbon supports a half-integer quanti- 
zation of the supercurrent to (n + l/2)4eAo/H [53] with 
a current-phase relation identical to the ordinary SQPC. 

The paper is organized as follows. In next section 
we introduce the model of the graphene SQPC and find 
the solution of the DBdG equation. In Sec. 3 we describe 
the boundary conditions for different edges and calculate 
the corresponding Josephson current. Sec. 4 is devoted 
to analysis of the obtained results. Finally in Sec. 5 our 
conclusion comes. 



2 Model and basic equations 

The geometry of the studied Josephson junction is shown 
schematically in figure [T] A ballistic graphene strip (N) 
connects two wide superconducting regions (S) which are 
produced by depositing of superconducting electrodes on 
top of the graphene sheets. Recent experiments have ver- 
ified the reality of this model [T81H9] . We assume the in- 
terfaces between the superconductors and the graphene 
strip are free of defects and impurities. The supercon- 
ducting parts are supposed to be heavily doped such that 
the Fermi wavelength Xf s inside them is very smaller 
than the superconducting coherence length £ and also 
the Fermi wavelength in the normal graphene strip Xp n- 
By the first condition mean field theory of supercon- 
ductivity will be justified and by second we can ne- 
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gleet the spatial variation of A inside the superconduc- 
tors close to the normal-superconductor (NS) interfaces. 
Thus the superconducting order parameter in the left 
and right superconductors has the constant values A = 
Aq exp(±i0/2), respectively, and vanishes identically in 
N. 

The low-lying electronic states in graphene can be de- 
scribed by the Dirac model as linear combinations of four 
zero-energy Bloch functions with slow varying envelopes 
(ipl, V>2, ip'i, Here ipi and ip2 are the wave function 
components on the two inequivalent atomic sites in the 
honeycomb lattice usually referred as pseudospin com- 
ponents. The primed and unprimed functions denote the 
two inequivalent Dirac points K and K' in the band 
structure of graphene known as the two valleys. The evo- 
lution of the envelope functions is governed by the two 
dimensional Dirac Hamiltonian [Zl[8l[9]: 

H = —i hv F (T Q ® a x d x + t z (gi g v d y ). (1) 

The Pauli matrices <ii and t\ act on the sublattice and 
valley degrees of freedom, respectively (with cto and To 
representing the 2x2 unit matrices). The eigenfunctions 
of the Dirac hamiltonian are the plane waves: 

ipK = e lkx+lqy (l,e la ,0,0), 
iP K , = e lkx+lqy (0,0,1, e~ la ), (2) 

in which the angle a = arctan^/fc) indicates the prop- 
agation direction. By treating the pseudospin similar to 
the real spin, one can easily see the pseudospin state 
(l,e ±m ) aligns parallel with the propagation direction. 
This allows one to introduce chirality, that is formally 
a projection of the pseudospin on the direction of mo- 
tion. The electron and hole have positive and negative 
chirality, respectively [10]. In ballistic systems chirality 
of states does not change and this conservation of chiral- 
ity leads to important effects in scattering phenomena, 
which Klein paradox is one of them pi?l[T5] . 

The superconducting correlations between the electron- 
hole excitations are described by the Dirac-Bogoliubov- 
de Gennes (DBdG) equation for the electron and hole 
four-component wave functions u and v |15j: 

(*■/■" „.*»-•-.) (:)=»■ < 3 > 

Here A is the superconducting pair potential matrix and 
is the time reversal operator. The excitation energy e 
is measured relative to the Fermi energy fj,. Each of the 
four blocks in Eq. [3] represents a 4 x 4 matrix, acting on 
two sublattices and two valleys spaces. The matrix forms 
of the pair potential, and the time reversal operators are: 

A = AT x ®a , (4) 

e = T X ® a z c. (5) 

in which C is the operator of complex conjugation. 



The singlet superconducting pairing always occurs 
between the time-reversed states. So in the case of the 
Dirac fermions in graphene, which the time-reversed states 
belong to the different valleys, the superconducting pair- 
ing is between the particles with different valleys as well 
as opposite spins and momenta. This cause the time re- 
versal symmetry to play a special role in the Josephson 
effect compared to that in the phase coherent transport 
in the normal state. Indeed the presence of a supercon- 
ducting electrode provides an intrinsic mechanism that 
couples phase coherently the electronic states belonging 
to the opposite valleys. So the dynamics of electrons is 
described by the full two- valley Hamiltonian of graphene 
as indicated in the DBdG equation. 

The phase difference <fi between the order parameters 
drives a supercurrent through the graphene strip which 
constitutes a weak link between two superconductors. 
This Josephson supercurrent is carried by the so called 
Andreev states which are formed in N region due to the 
successive conversions of the electron-hole excitations to 
each other (Andreev reflection) at the NS boundaries. 
For a short junction of L <C ^ and at zero tempera- 
ture the Andreev bound (discrete) states with energies 
|e| < Ag have the main contribution to the Josephson 
supercurrent [24]. In this case we can neglect the con- 
tribution from the continuous states above the super- 
conducting gap. We obtain the energies e(4>) of Andreev 
bound states by solving the DBdG equations with appro- 
priate boundary conditions which are described below. 
The Josephson current is then obtained from the formula 



between the Josephson current and the derivative of the 
free energy F with respect to the phase difference. 

Inside N the solutions of the DBdG equation are inde- 
pendent electron and hole-like wave functions which are 
classified by the 2-dimensional wave vector k = (fc, q) 
with the energy- momentum relation e — ftv^|k|. For a 
finite width W the transverse momentum q is quantized 
by imposing the boundary conditions at the edges. These 
transverse boundary conditions are different for the dif- 
ferent edge types which a graphene strip can have. 

The solutions inside the superconductors are rather 
mixed electron-hole excitations, called Dirac-Bogoliubov 
quasiparticles. The interfaces scatter the particles be- 
tween the two neighboring regions. However, the exci- 
tations inside N with energy inside the gap cannot en- 
ter the superconductors and two kinds of processes, An- 
dreev and normal reflection can occur for them. In gen- 
eral to match the solutions inside the different regions, 
one needs the scattering matrices of the interfaces. But 
assuming the ideal NS boundaries electron-hole scatter- 
ing can be described by a longitudinal boundary condi- 
tion between the electron and hole wave functions which 
reads [Tl] . 

u = t x <g> e -^+ l ' 3n -' T v , 13 = arccos(e/Z\ ) (7) 
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where <p is the phase of A in S and n is the unit vector 
perpendicular to the interface pointing from N to S. 

Assuming elastic processes, scattering at the inter- 
faces does not change the transverse momentum q and 
there is not any scattering between the different modes. 
So the Andreev states will be linear combinations of the 
wave functions from the same mode n with certain trans- 
verse momentum q n . For each type of the edges by im- 
posing the longitudinal boundary conditions at the in- 
terfaces to a general linear combination of the electron 
and hole states of mode n we reach to an equation for 
the Andreev energies e n (</>). Throughout the paper we 
will consider the limit Aq <C [i where the retro- Andreev 
reflection dominates [15] . In this regime the longitudinal 
momenta of the electron and hole k n (±e) can be approx- 
imated with k n (0) independent of the excitation energy 

H3|. 

3 Josephson current of the nanoribbons with 
different edges 

3.1 Smooth edges 

A graphene strip with smooth edges corresponds to an 
infinite mass confinement which is described by the fol- 
lowing boundary conditions 30 , 

V>i(0) = -^a(O), MW) = MW). (8) 

These conditions do not mix the two valleys and we can 
treat them independently. By this conditions the trans- 
verse momentum is obtained to have the quantized val- 
ues q n = (n + t})tt/W. We note that this quantization 
condition for the Dirac particles confined by an infinite 
mass differs from the corresponding one for the ordinary 
normal electrons by an offset of 1/2. This originates from 
the Berry phase 4>b — "R m monolayer graphene which 
cause 7r phase shift in the boundary conditions as well [3[ 
I30j . Considering a linear combination of the plane wave 
solutions (Eq. [2]) with opposite transverse momenta and 
applying the boundary conditions we can find the modes 
inside the strip. The electron modes for each valley are 
indicated by the wave functions of the form: 

u k = e lkx [sm(qy - a/2) , sm(qy + a/2)]. (9) 

The hole wave functions v k for each particular mode are 
the time reversed of the corresponding electron mode. 
Andreev bound states for each mode will be a linear 
combination of u± k and v±k- 

fu A \fau k + bu_ k \ 
\v A J \a'v k + b'v- k J 

To form an Andreev state this general form should sat- 
isfy the longitudinal boundary conditions (Eq. [7]) at the 
NS interfaces (x = 0,L). These conditions result in four 
independent homogeneous equations for four coefficients. 



By solving these equations we obtain the quantized An- 
dreev energy as 

en(4>) = A y/l-t n sin 2 (0/2), (11) 

with t n = [1 + (q n /k n ) 2 sin 2 (fc„L)] _1 . Comparison with 
Ref. [13] shows that t n is indeed the transmission prob- 
ability for a ballistic smooth-edge strip of graphene be- 
tween two heavily doped electrodes in the normal state. 
This type of relation between the Andreev energies and 
normal-state transmission probabilities t n is a universal 
form in the short junction limit (£ ^> L) [33]. By means 
of this method we can also find the transparency of each 
mode and then calculate the normal state conductance 
(Gn) using the Landauer-Buttiker formula at zero tem- 
perature: 

n 

The factor 4 at above relation accounts for the two-fold 
spin and valley degeneracies. Now from the Andreev en- 
ergies and also Eq. [S]the Josephson current at zero tem- 
perature can be found. 

Let us concentrate on the transmission coefficients 
and their dependence on the size and geometry of the 
junction. We note that in general both types of oscil- 
latory (real k n ) and evanescent (imaginary k n ) modes 
participate in the transport. The number of propagating 
modes is limited by the width as compared to the Fermi 
wavelength. Indeed only the modes with n < Nyy = 
fiW/nhvp — 1/2 can propagate and all the upper states 
are evanescent. Such a distinction between the modes 
exists in an ordinary QPC but with an important differ- 
ence. In the QPC the transparencies depend on the lon- 
gitudinal kinetic energy En — h 2 k 2 l /2m as t„ = 1/(1 + 
e - £ Jl/C^ j n w hich £ is a function of the geometry and 
size of the channel |31]f32] . Because of the quadratic form 
of the dispersion relation and high Fermi energies in 
comparison with £, we have |JSj[/C| 3> 1 for all modes. 
Consequently all the propagating modes are completely 
transmitted and the evanescent modes have very small 
transparencies. But since in graphene very small Fermi 
energies are accessible and the energy depends linearly 
on the momentum the evanescent modes can have a sig- 
nificant contribution in the transport. In particular near 
the Dirac point fiL/fwF *C 1, the number of propagating 
modes are very small and the evanescent modes have the 
main contribution in the transport. This is the unique 
property of the graphene Josephson junctions in which 
very small Fermi energies are accessible. In this limit the 
transmission coefficients through ballistic graphene are 
approximately, 

tn = cosh~ 2 (g„L), (13) 

which have the same form as transmission eigenvalues 
of the diffusive transport through a disordered contact 
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3.2 Armchair edges 

From Fig. Q] (lower panel) one may see that at the arm- 
chair edges the termination consists of a line of dimers. 
Therefore the wave function amplitude vanishes on both 
sublattices at the two edges. This leads to the boundary 
conditions, 

at the both edges. In contrast to the smooth edges, arm- 
chair edges mix the two valleys and causes scattering 
of a wave propagating on a certain valley to the oppo- 
site valley. The boundary conditions leads to the quan- 
tized transverse momenta q n = nit/W similar to the 
usual confinement of normal electrons. Due to the val- 
ley changing property of these edges, the electron wave 
function of each mode inside the nanoribbon is a mixture 
of the two valleys with the form 

u = e lkx [e lqv {\, e ia ), e~ lqy {l,e la )]. (15) 

This wave function has the same chirality for both val- 
leys contributions, as a result of the symmetric bound- 
ary conditions. Also one should note that in contrast to 
the smooth edges, the wave functions u(q) and u(— q) 
represent distinct solutions in which the components of 
the two valleys are interchanged. As a result the low- 
est mode n = has not the two fold valley degeneracy 
of the higher n = 1,2, ... modes. This zero mode with 
zero transverse momentum can exist inside the channel 
even if its widths becomes smaller than the Fermi wave- 
length. Like the case of the smooth edges the Andreev 
state wave function corresponding to a given mode n is 
the combination of u±k and v±k- But there is a difference 
with the smooth case. This time ua and va have four 
nonzero components. Applying the longitudinal bound- 
ary conditions leads to 8 equations for the coefficients 
a,a',b,b' in Eq. 1101 However because of the symmetry 
between the contributions from the two valleys, only 4 
of them are independent. We find that the Andreev en- 
ergies have the same form as in the smooth strip (Eq. 
ITT]) . So we would expect to reach similar results for I c 
and Ga? as well. 

3.3 Zigzag edges 

We now consider the case of the zigzag edges. It is clear 
from Fig. [TJ (upper panel) that atoms at zigzag edges 
belong to different sublattices. The boundary conditions 
are = ^ = at y = and ip 2 = = at V = W 
[34] . which lead to the transcendental relation, 

sm(qW)/qW = ±hv F /{v + e)W, (16) 

for the transverse momentum q. In contrast to the previ- 
ous cases this relation couples q to the excitation energy 
e. Of course in the regime of our interest (|e| < Aq -C fJ>) 
we can neglect e at Eq. 1161 This equation has a finite 



number of solutions depending on the value of [iW/hvp. 
Indeed all possible values of the transverse momenta are 
smaller than fj,/Hvp and thus the longitudinal momenta 
always have real values. This means that all the modes 
inside the nanoribbon with zigzag edges are propagat- 
ing. For /jW/hvF < 1 there is an imaginary solution q 
labeling an evanescent mode in the y-direction. In the 
interval 1 < [iW/hvp < 3ir/2, the transcendental rela- 
tion has a single oscillatory solution. For larger fiW/fiVF 
the number of the solutions increases by two whenever 
the width is increased by a half of the Fermi wavelength 
Xf n = hvp/lJ" We classify the wave functions to two 
groups according to ± sings in the transcendental rela- 
tion. To find the wave functions we consider the general 
combination of the plane wave solutions: 

4, = A4, q K + B^ K q + A'4, q K>+B'^- K q . (17) 

The boundary condition at y = leads to A + B = 
A' + B' = 0. So the wave function must be as below: 

(Asm(qy) \ 
AMqy + a) ( } 

A'sm(qy) I v ; 

A' sin(<?y — a)/ 

Now to satisfy the second condition at y — W there 
are two ways. We can have either s\a.(qW + a) = and 
A' = or sin(qW — a) = and A = 0. We consider the 
solutions belonging to the group with sm(qW)/qW = 
+Hvf/ fiW. If k is positive the first way for satisfying the 
boundary conditions will be possible and if k is negative 
then the second way is acceptable. So the electronic wave 
functions of the first group (+ sign) have the forms, 

u + = e ikx \s\n{qy) , sm{qy + a) , 0, 0] , 
U~ = e' ikx [0,0 1 sm{qy) 1 sm(qy ~ a)}, (19) 

which, respectively, describe a right-going wave on the 
valley K and a left-going wave in the other valley K' . It 
can be easily seen that the wave functions of the second 
group (— sign) are obtained from the above functions by 
the replacement k —> —k. We note that for each mode the 
zigzag nanoribbon operates as a valley filter for the waves 
with mixed valley components [23j . Since the valley in- 
dex and q are conserved upon a normal reflection, the 
filtering property prevents the electrons from such scat- 
tering. Consequently the wave function of each bound 
state in N will be a combination of and v~( + \ 

From the longitudinal condition (Eq. [7]) we obtain that 
the energy of the Andreev state does not depend on the 
transverse mode n and is given by, 

e = Z\ o cos(0/2), (20) 

which is the same as the Andreev energies of an ordinary 

SQPC gl]. 
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1 2 3 4 5 6 
W/L 

Fig. 2 The critical current (a) and the product I c Rn (b) 
dependence on the relative width of junction at a typical 
small value fj,L/hvF = 0.1 for smooth and armchair strips. 



4 Results and discussion 

Let us first consider the behavior of the Josephson cur- 
rent near the Dirac point when [iL/Tinf *C 1. Fig. [2] 
shows the critical current and the product I c Rn vari- 
ations with W/L for the smooth and armchair edges 
at a typical small Fermi energy fiL/hvp = 0.1. We see 
I c has a linear dependence on W/L at large W/L fl2k) 
which is the result of diffusive-like transport in the bal- 
listic graphene [2]. By decreasing W/L, the transmis- 
sion probabilities of the modes through N decrease and 
the critical current shows a monotonic decrease with- 
out any quantization. For a narrow strip W < L the 
evanescent modes transparencies are vanishingly small 
and I c reduces to a constant minimum value. For the 
smooth edges there is no propagating mode when W < 
{tt / 2)hv p / ji and the minimum supercurrent is vanish- 
ingly small. However in the case of the armchair edges 
a lowest nondegenerate mode with q = always can 
propagate irrespective of the width of the junction. This 
zero mode results in a bound state carrying a nonzero 
residual supercurrent I c = eA /H for W < L. (Fig. Oi) 

For a wide contact the normal state conductance Gn 
has the same linear dependence on W/L as I c and thus 
the product I c Rn reaches a constant value 2.08Z\ /e in 
the limit W S> L(see Fig. (2b) for both types of the edges. 
Decreasing the width the difference between the Joseph- 
son coupling constants of the junctions with two differ- 
ent edges appears. For smooth (armchair) edges I c Rn 
decreases (increases) monotonically with lowering W/L 
and reaches to a minimum (maximum) of (l/2)-7n4o/ e 
(irAo/e) for a narrow junction of W < L (see Fig. |2Jd). 

As indicated in Fig. [3] the current-phase relation of 
finite- width junction differs significantly for the armchair 
and smooth cases specially for very narrow strips. How- 
ever at the limit of very wide junctions W L for both 



W/L=100 

W/L=1 

W/L=0.01 




(b) 





W/L=100 




W/L=1 




— ■ W/L=0.01 


Armchair 




- 





71 

3> 



2tz 



Fig. 3 The current-phase relation for junctions with smooth 
and armchair edges at different W/L values for fj,L/hvF = 
0.1. 



types, the current-phase relation is like a diffusive metal 
junction fl6 | fl7 } : 



eA 2W 
H ttL 



cos(</>/2)arctanh(sin(</>/2)). (21) 



On the other hand a very narrow strip W <C L with 
smooth edges has a sinusoidal current-phase relation /(</>) 
I c sin((/)) similar to a tunnel Josephson junction. This is 
expected as in this limit the transmission probabilities 
are very small and the graphene ribbon behaves as a tun- 
nel junction. For a narrow armchair ribbon the current- 
phase relation reads 



(22) 



which is the characteristic relation for the SQPC. This is 
also understandable because even in a very narrow arm- 
chair strip there is a propagating mode with complete 
transmission which makes it similar to the SQPC. The 
current-phase relation of smooth and armchair junctions 
are shown in Fig. [3] for different W/L. 

Now let us consider the Josephson current at a higher 
Fermi energy far from the Dirac point. Fig. 2^ and Ujo 
show the dependence of I c and I c Rn on the junction 
width, respectively, for (iL/hvp = 5. In this regime both 
the propagating and evanescent modes contribute to the 
supercurrent. Roughly speaking the amplitude of the 
contribution from the evanescent modes is proportional 
to W/L, however, the propagating modes contribution 
is proportional to their number Nyy ~ [iW/itHvf- When 
/iL/hvp <C 1, the evanescent modes play the main role 
in the transport and thus the transport is diffusive-like 
and both the critical current and conductance are pro- 
portional inversely to length of the junction L. On the 
other hand, when ^L/hvp ^ 1 the propagating modes 
have the major contribution. This cause that the enve- 
lope of the critical current depends linearly on /j,W/hv F . 
In addition to this overall increase with jiW/hvp the 
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4 8 12 16 20 

(iW/hv F 

Fig. 4 The critical current (a) and the product I c Rn (b) 
dependence on the relative width of junction at a typical 
large value [iL/hvF = 5 for smooth and armchair strips. 



critical current undergoes a series of peaks which become 
smoother by increasing fiW/hvF- Each peak (jump) sig- 
nals addition of a new propagating mode in the transport 
when the width increases by a Fermi wavelength. Both 
the critical current and the product I c Rn have an oscil- 
latory behavior versus [iW/Hvf with a period of 2n (Fig. 
[4j . We should note that due to a nonzero transverse mo- 
mentum of the lowest mode (n = 0) in the junction with 
smooth edges (while the zero mode of the armchair-edge 
nanoribbon has zero transverse momentum), the oscilla- 
tions are shifted by half a period with respect to those 
of the armchair edges case. Also comparing the results 
of Fig. [2] and [H we see that the limiting values of I c and 
I c Rn at W < L is always the same and does not de- 
pends on the value of the chemical potential /i. Therefore 
we find that for these two types of the edges there is no 
sharp quantization and stepwise variation with width for 
I c and Gn- This is again a unique property of graphene 
with smooth and armchair edges in which the evanescent 
modes can contribute significantly in the transport. 

Finally let us analyze the Josephson supercurrent for 
the zigzag ribbon. As we argued in section [3] in this case 
the valley filtering nature of the wave functions in N pre- 
vents the normal reflections at NS interfaces. In addition, 
because of the transcendental relation of the transverse 
momenta which confines their values by the Fermi en- 
ergy /z, there is no evanescent mode with imaginary k 
and the situation is similar to an ordinary SQPC. The 
critical supercurrent shows a step-wise variation with 
fiW/hvp but with the following important anomalies 
(see Fig. [5]). In contrast to an ordinary SQPC [2"4"I|2"5] . 
the width of the first step (3/4) is bigger than that of the 
higher steps (1/2). The extra width is the contribution 
of the single y-directional evanescent mode (imaginary 
q) for W < Hvp/fj,. Also the height of the first step is 
1 /2 of the height of the higher steps which itself is four 




o 1 

HW/hv F 

Fig. 5 The critical current vs. relative width of junction 
which has zigzag edges. 



times bigger than cAq/H the supercurrent quantum in 
an ordinary SQPC. Therefore the supercurrent through 
a zigzag graphene nanoribbon is half-integer quantized 
to (n + l/2)4eAo/H. The effect resembles the conduc- 
tance quantization in graphene strips with zigzag edges 
[T2] and also the half-integer Hall effect [HOSE] m mono- 
layer structures. From expression of the Andreev ener- 
gies for the zigzag junction (Eq. |20|) we find that the 
current-phase relation is given by Eq. [2U similar to an 
ordinary SQPC with length smaller than £. 



5 Conclusion 

In conclusion we have investigated the Josephson effect 
in a short graphene nanoribbons of width W and length 
L -C £o connecting two heavily doped superconducting 
electrodes. Within a Dirac-Bogoliubov-de Gennes for- 
malism, we have found that the variation of the Joseph- 
son supercurrent versus the width W is drastically dif- 
ferent for different types of edges. In the smooth and 
armchair nanoribbons with low concentration of the car- 
riers the critical supercurrent I c decreases monotoni- 
cally by decreasing W/L. For a narrow strip W < L 
with armchair and smooth edges, I c takes constant min- 
imums eA^/H and 0, respectively. The Josephson cou- 
pling strength given by the product I c Rn has been found 
to have a minimum (l/2)7rZ\o/e (maximum irAo/e) for 
a narrow junction W < L with smooth (armchair) edges 
and increases (decreases) monotonically with W/L to 
reach the wide junction value 2.08 Aq/c We have also 
seen the different dependence of the current-phase rela- 
tion on W/L for junctions with the two different edge 
types. 

For the higher concentration of the carriers, the over- 
all monotonic dependence of I c acquires a series of peaks 
with distances inversely proportional to the chemical po- 
tential [i. Correspondingly the product I c Rn shows an 
oscillatory behavior with fjW/hvp with a period of 2-7T. 
The oscillations for smooth and armchair edges are phase 
shifted by tt. 
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For the zigzag edges the results are quite different 
due to the special quantization relation of transverse 
momenta and the valley filtering property of the electron 
wave functions. We have found a step- wise variation of I c 
versus fiW/hvF, implying a half-integer quantization of 
the supercurrent to (n+ l/2)4e Aq/H. The current-phase 
relation is found to be similar to an ordinary supercon- 
ducting quantum point contact. 
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